Predicted signatures of p-wave superfluid phases and Majorana zero modes of 

fermionic atoms in RF absorption 
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We study the superfluid phases of quasi-2D atomic Fermi gases interacting via a p-wave Feshbach 
resonance. We calculate the absorption spectra of these phases under a hyperfine transition, for both 
non-rotating and rotating superfluids. We show that one can identify the different phases of the 
p-wave superfluid from the absorption spectrum. The absorption spectrum shows clear signatures of 
the existence of Majorana zero modes at the cores of vortices of the weakly-pairing p x + ip y phase. 
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The observation of p-wave Feshbach resonances in cold 
fermion gases [l| opens up the possibility of the exper- 
imental realisation of unconventional superfluid states 
with non-zero pairing angular momentum^, 0, 0|- It 
has recently been predicted [H, H| that a one-component 
fermion gas interacting via a p-wave Feshbach resonance 
shows a series of p-wave superfluid phases, including &p x 
phase, and both strong- and weak- pairing phases of 
p x ± ip y symmetry. The strong- and weak- pairing phases 
differ in the nature of their excitation spectra, and are 
separated by a topological phase transition [g, LZt] - 

Among these, the weak-pairing p x ± ip y phase has re- 
cently received considerable theoretical attention. The 
vortex excitations of this phase are predicted to hold 
gapless Majorana fermions on their cores, which should 
give rise to non-abelian exchange statistics in a quasi 2D 
geometry 0, Such excitations may have their uses in 
topologically-protected quantum computation 0, @, E|. 
However, to our knowledge, there as yet are no exper- 
imental observations of these very interesting physical 
properties in any realisation of a p x ± ip y superfluid. De- 
veloping methods to experimentally probe the properties 
of p-wave superfluid phases in cold atom systems is then 
of high and timely interest. 

In this Letter we show that features in the RF absorp- 
tion spectrum of a fermionic gas provide unambiguous 
signatures of the different p-wave superfluid phases. We 
show that the weak-pairing phases can be identified from 
the strong-pairing phases in the absorption spectrum of 
the superfluid at rest. Most significantly, we focus on the 
Px + ipy phase and show that for a rotating fluid at weak- 
pairing, the absorption spectrum has clear and striking 
signatures of the zero energy Majorana fermion modes 
on the cores of the vortices. This spectrum is therefore 
a direct probe of the physics underlying the proposed 
non-abelian exchange statistics of this phase. (These ex- 
perimental signatures do not require manipulation of po- 
sitions of vortices or the use of other localised probes [9(.) 
Furthermore, we show that the absorption spectrum for 
the rotated p x + ip y weak-pairing phase depends on the 
sense of the rotation. This dependence is a direct in- 
dication of the time-reversal symmetry breaking in this 



phase. 

We consider a gas of identical fermionic atoms in an 
internal state that we denote ||), interacting via a p- 
wave Feshbach resonance. The phase diagram of the 
system has the same qualitative form in wide Q and 
narrow Q resonance limits. Depending on temperature, 
and on the detuning and anisotropic splitting of the Fes- 
hbach resonance, there appear superfluid phases with p x 
and p x + ip y symmetries. For p x + ip y , there is a transi- 
tion between the weak-pairing phase (chemical potential 
fi > 0) and the strong-pairing phase (/x < 0) as the res- 
onance is swept from positive to negative detuning^. 
We propose to probe the properties of these phases by 
measuring the absorption spectrum, under excitation of 
an atom in state ||) to a new internal state |f) via the 
perturbation 



pert 



n 



d 2 re lqr ~ 



l ijA{r)%j;i(r) + h.c. 



(1) 



where ij sets the Rabi frequency of the transition. We 
consider \\) — > |f) to be a hyperfine transition, driven 
by RF absorption [loi [llj or a two-photon transition 
measurement of the population of excited (or remain- 
ing) atoms is achieved by separate absorption on an elec- 
tronic transition. We denote the energy splitting between 
the states ||) and |f) by E g ; for a hyperfine transition, 
Eg r~j 10 9 Hz, and we set the net photon momentum trans- 
fer to zero, q — > 0. The "["-atoms will not participate in the 
resonant p-wave interaction with the |-atoms, so for the 
most part we shall consider the T-atoms to be free. How- 
ever, we shall discuss the effects of s-wave interspecies 
interactions. We shall focus on a uniform gas in a quasi- 
two-dimensional geometry, that is with confinement ap- 
plied in the z-direction such that the associated confine- 
ment energy is large compared to the Fermi energy. This 
simplifies the analysis, allowing analytical calculations of 
the absorption features. The three dimensional case is 
left for a future work. 

When the interaction between atoms may be neglected 
(far from the Feshbach resonance) and they form a sim- 
ple Fermi gas, a transition of an atom from (J.) to ||) re- 
quires an energy of E g , leading to a delta- function peak 
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of the absorption at u> = E g . (We choose units for which 
h = 1.) In contrast, we find that in the superffuid phases 
(close to the resonance) the absorption differs strongly 
for weak- and strong-paired superfluids. In the weak- 
pairing phase the energy E g becomes a threshold energy 
for absorption, with non-zero weight at uj = E g and a 
continuous spectrum above this (see Eq. ((4])). In the 
strong pairing phase the threshold is shifted to E g +2|/x|, 
with weight that is zero at E g + 2\fi\ and grows linearly 
with energy. The discontinuity in the absorption spec- 
trum in the weak-pairing phase unambiguously identifies 
it from strong pairing. 

In the weak-pairing p x ± ip y phase, the introduction 
of vortices by rotation leads to the creation of a set of 
gapless modes on the vortices. In the presence of a vor- 
tex lattice we find that a series of equally spaced sub- 
gap peaks will be present in the absorption spectrum, see 
Fig. pQ. These peaks are a result of transitions from 
the zero energy core states of |-atoms to various states 
of f-atoms, and have a weight that is linearly propor- 
tional to the number of vortices Ny- The peaks start 
at ~ E g — fi, and their weight decays monotonically in 
energy, becoming small well before the onset of the con- 
tinuum at E g — fi + A (A is the superfluid gap). For 
a dense vortex lattice, in which there is tunneling be- 
tween the Majorana modes, the peaks broaden to reveal 
a lineshape showing a van-Hove singularity, see inset of 
Fig. [Tj. For the same weak-pairing phase, with vor- 
tices of the opposite sense of rotation, we find that the 
absorption in the Majorana modes again gives rise to a 
set of narrow absorption peaks. The peaks still start at 
~ E g — fi, but now their weight rises linearly in energy, 
passing through a maximum well before the continuum. 
As we shall discuss below, the difference in spectra for 
the two senses of rotation is a signature of the fact that 
Px + ipy breaks time-reversal symmetry. 

We now provide the details of the calculations. The 
rate of excitations between the two hyperfine states is 
governed by the Golden Rule 

r[w] = 2^M 2 Y, \M ab \ 2 S(E Ua + E u - u), (2) 

a b 

where E^ a (-&]■, ft) is the energy required to produce an 
|-excitation {] -excitation), and M a b is the dimensionless 
matrix element (measured in units of 77) for producing 
excitations with quantum numbers a and b respectively. 
In the following we shall calculate T[oj], first for a system 
at rest, then for a rotating system. 

We describe the superfluid phases within the 
Bogoliubov-de-Gennes (BdG) approach. For a uniform 
system, the Hamiltonian is given by 



H — i?0,J. + -ffpaii',1 + #0,T 
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where H Q>S = j d 2 k^ s (k) - fij tfj s (k), 
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FIG. 1: Peak distribution associated with the Majorana 
modes on the vortices of the weak-pairing p x + ip y phase. The 
first peak will be at E ~ E g — p., and other peaks will follow 
at spacings lo c apart. Continuous absorption will be mea- 
sured starting at E ~ E g . For the opposite sense of rotation, 
the weight of the peaks starts linearly and passes through a 
maximum. The inset shows the shape of a single peak for a 
square (dashed line) and triangular (solid line) lattice, show- 
ing a van-Hove singularity at E — 2t (P is measured in units 
of 2n\ V \ 2 NvOo/t). 



# P air,i = / d 2 k^l(k)v A (k x + ik y )^[(-k) + h.c. The 
values of fx and u A may be obtained from the mean-field 
calculations of Refsy, |j. The Hamiltonian for .[-atoms 
is diagonal in the Fock space of au — Wfc"0i.fc + v k^l 



and cr k = u^tp] k + u^-0l,-fc where Uk,Vf, are solutions 
of the BdG equation, |w fe | 2 = \ \\ + - fij /E iM 

and E^k = 



and 



\Vk\ 



v\k 2 0,0]. Treating the excited t-atoms 
as free, E^k = fc 2 /2m — p, a summation over k leads to 
the absorption spectrum [I9( 
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9(5u-2\fi\), n<0, 
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where <S* is the area of the sample and 5u> = uj — E g . Non- 
zero interactions of the excited "["-atom with the J,-atoms 
lead to a (mean-field) shift of both of these spectra • 

We now turn to discuss the system under rotation, in 
a regime where a large vortex lattice is formed, Ny ^S> 1. 
(The rotation frequency is then close to the trapping 
frequency lot-) We first examine a superfluid in the weak 
pairing p x + ip y phase, and concentrate on the regime 
f2 < fi < mv\. 

Projected to the low-energy sector of the Majorana 
modes, the Hamiltonian for the |-atoms becomes a tight- 
binding Hamiltonian of the form[l2T| 

H i = H 0,l + #pair,4, = it Y Simile > ( 5 ) 

(ij) 

where = ±, in a way that corresponds to half a flux 
quantum per plaquette for a square lattice, and a quar- 
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ter of a flux quantum per plaquette for a triangular lat- 
tice. Here 7$ represents a Majorana fermion localized at 
Ri. The parameter t describes the tunneling between 
the Majorana modes of nearby vortices; this tunneling is 
assumed small, t -C ji, which is valid if the separation 
between vortices, I = ^/l/(2mO), is large compared to 
the spatial extent of the Majorana mode, tq. The Majo- 
rana wave function g{r) has an oscillating contribution 
within the vortex core, and a decaying part outside the 
core. When /1 <C mv\, the size of the core can be es- 
timated as r coro <~ log ' — , with vf being the 
Fermi velocity of the underlying fermionic gas, while the 
decay length outside the core is tq — v & /p. In the limit 
of small [i we have tq 3> r COIC , such that we may approx- 
imate g(r) « "^73==! with r being the distance from the 
center of the vortex 11311 . The analysis of the Hamiltonian 
([5]) is presented in |l2l . [Til ] . For a square and triangular 

lattice the spectra are E n ,k = 2t^/ 'sin 2 (ak x ) + sm 2 (ak y ) 

and E & j- = y/2t\J 3 + cos(2afc x ) — 2 cos(ak x ) cos(V3ak y ) 
respectively (l2^. 



The spectrum of the "[-atoms in the rotating frame 
depends on the strength of the s-wave interspecies in- 
teractions. Since this depends sensitively on the atomic 
species used, we shall consider two limiting cases: For 
vanishing interspecies interactions, the "("-atoms arrange 
into Landau levels due to rotation, with the cyclotron 
frequency lo c — 20, ~ 2n x 10 2 Hz. (The confinement 
frequency for "[- and |-atoms is assumed the same.) For 
strong interspecies repulsion, the low energy states of the 
"["-atoms are localized at the cores of vortices where the 
density of J,-atoms is low, and can be described by a 
tight-binding Hamiltonian. The rotation then translates 
to a magnetic flux threading the lattice plaquettes (the 
Azbel-Hofstadter problem [15j|). 

Starting with the case of zero interspecies interactions 
(Vf = 0), the Hamiltonian i?o,t can be diagonalized by 
expanding Vf( r ) = Y, np ^n P {r)c np , where (j> np {r) are 
wavefunctions for the n'th Landau level, and p is some 
quantum number related the degeneracy of the Landau 
level. The rate of energy absorption is determined by 



v y 



where gi(r) — g(r — Ri). The functions A"^ k are listed 
in Refflil. z(i) = 1, . . . ,v numbers the vortices contained 
in a unit cell, and a = ±1, . . . ,±u/2 is a band index. 
Since p enumerates degenerate states, we can now sum 
over it using the addition theorem for Landau levels [l6| . 
We keep only the leading term Ri — Rj. Finally we 
expand the phases in Eq. ([6]) to first order in r — r', 



l + 2n 1 



Thus 



and obtain Y^ P \ M ak,n P \ = 
the integrated absorption into the n'th Landau level falls 
monotonically with increasing n. Plugging the matrix 
elements into Eq. ([2]), we arrive at the result 



r[w] = 2n l -^-N v J2e n F 



t 



t 



(7) 



where r 



E g -fi+uj c (n- 



ft) [l + 2n(?)' 
1/2), and F is a dimensionless function proportional to 
the density of states of the band formed by the Majorana 
operators, which is nonzero over a range of order unity. 
It is explicitly calculated below. 

Using the results above, assuming t <C u> c , we can find 
the integrated absorption over a single Landau level 



d(Lu-u n )r = Tr\n\ 2 N v G r 



(8) 



The integrated absorption over all frequencies for excita- 
tions from the zero modes can be found using the com- 
pleteness of the Landau level functions ^2 np \ M a k. np \ 2 = 
\, hence / 6ljY\u\ = ix\n\ 2 N v jl. 

When interspecies interaction is strong, the "[-atoms 
feel a periodic potential Vf due to the vortex lattice struc- 
ture. Due to depletion of |-atoms at a vortex core, the 
potential Vf possesses there a minimum, and may hold 
several bound states <fi n { r ) f° r T-atoms. Several tight- 
binding bands will be formed, labelled by n. The "["-atoms 
feel the same flux per plaquette as the Majorana fermions 
do. Therefore, the Hamiltonian for a single band of tight- 
binding atoms in the rotating frame may be written as 

H a A = KYl c Ul Cn ^ + it 'n^2 s ij c l,n c n^i ( 9 ) 
« (ij) 

where t' n is the tunneling matrix clement for the "[-atoms, 
and E' n is an on-site energy. Diagonalizing the Hamilto- 
nian, one finds that the spectrum of the "[-atoms has the 
same fc-dependence as that of the Majorana fermions. 
The resulting absorption is given by Eq. (0) with 0„ = 



t + t' 



I / d 2 rg(r)4> n {r)\ 2 , u n = E g - /z - 

We now proceed to calculate the function F appear- 
ing in Eq. Q. This function is F[E/t] = tp[E]/N v , 
with p(E) being the density of states of the Majorana 
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fermions. Starting with a square lattice, F[x] is 



Fn\x] 



Si- 



dy 



(2tt) 2 Jc \Jx 2 - y 2 a/4 - x 2 + y 2 \/4 - y 2 



.(10) 



For < x < 2, C 
[v^ 



= [Q,x], while for 2 < x < 2^2 
We note the diverging density of 



C = LV^ 2 -4,2j 
states at E = 2t, and the discontinuity at the top of 
the band E — 2\f2t. For a triangular lattice, we again 
find a logarithmic divergence near E = 2t and disconti- 
nuities in the response at the top and the bottom of the 
band, E = 2y/3t and E — y/3t respectively. 

For the opposite sense of rotation, the vortex lat- 
tice consists of an array of vortices of opposite circula- 
tion. These vortices also have zero modes [17j. However, 
the Majorana modes now have a non-trivial dependence 
around the vortex center, gt(r) oc e -larg (' Ri \ (We 
again assumed tq 3> r COIC .) As a consequence, the inte- 
grated weight of the absorption from the Majorana mode 
into the n'th Landau level is proportional, for small n, to 



s>r-£(?) 



(ii) 



The weight passes through a maximum at a frequency of 

2 

about -cS-t- above the onset at E a — a. 

rav ^ y r 

Thus we find that in the weak pairing p x + ip y phase, 
the spectrum depends strongly on the sense of rotation. 
This dependence may be probed by preparing the system 
in the p x + ip y phase and comparing spectra for rotations 
of the two senses. Alternatively, since what matters is 
the relative sense of rotation as compared to the internal 
angular momentum of the order parameter, the depen- 
dence could be tested by comparing the spectra of the 
Px + ipy and p x — ip y phases for the same sense of ro- 
tation. In a stationary fluid, the two phases are equally 
likely to be formed 2l|. If for a given sense of rotation 
absorption spectra of two types are randomly observed 
for different realizations of nominally the same parame- 
ters, it would be a direct demonstration of time-reversal 
symmetry breaking. We note that the results may de- 
pend on how the state was formed (apply rotation and 
then cool, or cool and then apply rotation) [22| . 

The observation of absorption that originates from Ma- 
jorana modes would also distinguish the time-reversal 
symmetric p x phase from the phase of p x + i(3p y , with 
a real (3 / 0, where time-reversal symmetry is bro- 
ken. The latter may appear for an anisotropic Feshbach 
resonance 0, The distinction arises, as we now show, 
from the presence of Majorana modes for all (3^0, and 
their absence for (3 — 0. Solutions to the BdG equations 
come in pairs of ±E. When a vortex carries a localized 
single zero energy state, as it does for /3 = 1, the only 
ways the energy of that state may shift due to a change 
of (3 are through mixing with other zero energy states (on 
other vortices, or on the sample's edge) or when the bulk 



gap closes. When vortices are far away from one another 
and from the edge, the zero energy state exists as long 
as the gap does not close, i.e., as long as j3 ^ [l8j][23j]. 

To conclude, our results show that the absorption spec- 
trum can be used to distinguish the various p-wave su- 
perfluid phases of an atomic Fermi gas, and to probe the 
breaking of time reversal symmetry. Furthermore, for a 
rotating superfluid in the weak pairing p x + ip y phase, 
the absorption shows a series of sub-gap peaks, of unique 
shape and distribution, which are the direct consequence 
of the zero-energy Majorana modes on the vortices. The 
most favourable parameters for the observation of these 
sub-gap peaks are [i <C mv\ (to avoid finite energy bound 
states in the vortex core), and temperature T <C \i (to 
avoid thermally excited quasi-particles which may wash 
out the absorption peaks). A practical way to observe 
the peaks would be by comparing measurements of the 
absorption for increasing rotation frequency. This will 
vary both the strength of the peaks (due to the increase 
of Ny) and their spacing (due to the increase of u c ). 
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